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Abstract 

We extend the resummation of dimensionally-regulated amplitudes to next-to-next-to-leading poles. 
This requires the calculation of two-loop anomalous dimension matrices for color mixing through soft 
gluon exchange. Remarkably, we find that they are proportional to the corresponding one-loop matrices. 
Using the color generator notation, we reproduce the two-loop single-pole quantities H^^^ introduced by 
Catani for quark and gluon elastic scattering. Our results also make possible threshold and a variety 
of other resummations at next-to-next-to leading logarithm. All of these considerations apply to 2 — > n 
processes with massless external lines. 
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I. INTRODUCTION 



The description of partonic hard scattering in quantum chromodynamics (QCD) is central 
to the analysis of final states at hadronic colliders. The calculation of cross sections for such 
processes requires a combination of virtual and real radiative corrections, organized according 
to underlying factorization theorems. This is the case for higher-order calculations to next- or 
next-to-next-to-leading order in (NLO, NNLO, . . . ). It holds as well as for resummed cross 
sections, in which selected corrections associated with soft and coUinear gluon radiation are 
organized, at leading, next-to-leading or next-to-next-to-leading logarithms (LL, NLL, NNLL, 
. . . ) to all orders in ag- 
in, both fixed-order and resummed calculations the coherence properties of soft gluon radiation 
play an essential role. An anomalous dimension matrix for inclusive wide-angle soft gluon radia- 
tion was introduced in Refs. , and computed to leading order for quark and gluon scattering 
processes in Ref. P]. The one-loop matrix of soft anomalous dimensions has been applied to 
the NLL threshold resummation of iet cross sections 

bined with resummed form factors, was shown in Ref. 
dimensionally- regularized partonic scattering amplitudes 'm e = 2 — D /2. In this paper we will 
show how to compute this matrix directly at two loops, from a relatively limited set of diagrams 
in the eikonal approximation, using Wilson lines, giving as an explicit example quark-antiquark 
scattering. 

The full analysis given below applies to any 2 n partonic amplitude in dimensional regu- 
larization. The two-loop soft anomalous dimension matrix allows the exponentiation of next-to- 
next-to-leading infrared poles, which appear in the combination a'^{l / £)^~^ in the exponent, a 
level equivalent to next-to-next-to-leading logarithms. The resulting resummed amplitudes can 
be expanded out to the two-loop order, and the poles in e can be compared to explicit two-loop 
scattering amplitudes, for example the basic 2^2 scattering processes jlfl. 



y and of distributions of event-shape 
^ At two loops, the same matrix, com- 
to control the single infrared poles of 



terms presented in Ref. 



were expressed in terms of the color-space notation and the organization of two-loop singular 
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12l |. Those poles 



14| . (Related work at one loop was performed in Refs. 



15|, 
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We will 
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verify that the expansion of the resummed amphtudes to two loops matches precisely the full in- 
taed pole *uctu.e of the known two-loop scattering amplitudes, including the smgle poles in e. 

Remarkably, we will find, as reported in Ref. \Xf\, that the two- loop anomalous dimension matrix 
is related to the one-loop matrix by a constant, the same constant, appearing in the DGLAP 
splitting kernel, that relates the one and two-loop anomalous dimensions for the Sudakov form 
factor. (The analogous matrix appears in the electroweak Sudakov corrections to four-fermion 
scattering, and has been extracted at two loops from the QCD four-quark scattering amplitude 
in Ref. The simplicity of this result will facilitate the development of practical resummed 

cross sections with color exchange at NNLL. 

This paper is organized as follows. The next section reviews the collinear and infrared factor- 
ization of exclusive amplitudes. In that section, we provide a new explicit scale-setting choice for 
the soft function, which is necessary to define the scales of logarithms in the relevant anomalous 
dimensions. The third section describes the expansion of the jet functions to two loops. Here we 
describe a new "minimal" reorganization of the factorized amplitude, to facilitate the comparison 
to fixed-order calculations. In the fourth section, we describe in detail the one- and two-loop 
calculations necessary to determine the soft anomalous dimension matrix, for the specific case 
of quark-antiquark scattering. Here, we will employ the eikonal approximation, and the scale- 
setting choice for the soft function from Sec. 2. We show that diagrams attaching gluons to three 
different eikonal lines either vanish, or represent the exponentiation of the one-loop soft matrix. 
We close Sec. 4 by generalizing these calculations to arbitrary flavors for incoming partons and 
arbitrary flavors and numbers of outgoing partons. To do so, we present the color-mixing anoma- 
lous dimension matrix in the color-space notation of Ref. Il4|. Finally, in Sec. 5 we employ 
this notation, along with results of Sec. 4 for the soft anomalous dimensions and known two-loop 
elastic form factors for quarks and gluons, to give the explicit form of the two-loop single pole 
terms in e, for arbitrary 2 n partonic processes in QCD. We show that these pole terms agree 
with the single-pole "H^^)" terms found in NNLO 2 -> 2 calculations |ld, uJ, UJ, M whose 
poles have been organized according to the formalism of Ref. [l4|. Our results also agree with 
the proposal of Ref. j^l for the single poles for the case of 2 ^ n gluon processes, which was 
based on the consistency of collinear factorization of amplitudes. We provide an appendix with 



3 



explicit forms of Sudakov anomalous dimensions, and two appendices illustrating calculations of 
soft anomalous dimensions using eikonal methods. The final appendix details the computation 
of a particular commutator of color-space matrices, which is needed to compare our results with 
the explicit NNLO calculations. 



II. FACTORIZED AMPLITUDES IN DIMENSIONAL REGULARIZATION 

Our considerations apply to 2 — > n scattering processes, denoted as "f", 

f : '^l) + /2(P2, ^2) /3(P3, ^^3) + /4(P4, + ■ ■ ■ + /„+2(p„+2, ^^+2) • (2.1) 

The labels refer to the flavor of the participating partons, each of momenta {pi} and color 
{rj}. The amplitude for this process, Ai^^\ is a color tensor with indices associated with the 
external partons {rj} = {ri, r2, . . . }. It is convenient to express these amplitudes in a basis of C 
independent color tensors, {ci)^^..-^, so that P. Il^ 

= l-Mf) , (2.2) 

[fi 

where the ket may be thought of as a vector Ai^l with C elements in the space of color tensors 
cj. We will analyze these amplitudes at fixed momenta Pi for the participating partons, which 
we represent as 

p, = QA, f3-=0, (2.3) 

where the are four-velocities, and where Q is an overall momentum scale. For the purposes 
of this analysis, and to compare with existing NNLO calculations, we take all of the partons 
massless, as indicated. To be specific, we may take Pi ■ P2 = ^ for the incoming partons in 
Eq. ()2.H) . which implies = s/2, but this is not necessary. 

In dimensional regularization (D = 4 — 2e), on-shell amplitudes may be factorized into jet, 
soft and hard functions, which describe the dynamics of partons collinear with the external lines, 
soft exchanges between those partons, and the short- distance scattering process, respectively. 



This factorization follows from the general space-time structure of long-distance contributions to 



elastic processes j21 
presented long ago 



. A formal proof for the case n = 2 in QCD (quark-quark scattering) was 
The general form of the factorized amplitude is 

xi/f , (2.4) 



where /i is the renormalization scale. Jt^^ is the product of jet functions for each of the external 
partons, as above denoted collectively by [f], S''^^ is the soft function, and is the short-distance 
function. For example, when the process is 1 + 2 — > 3 + 4, the product of jet functions is 

^ ^ ^ j=l,2,3,4 ^ ^ ^ 

Construction of the soft and jet functions requires the specification of at least one independent 
momentum scale, Q', which plays the role of a factorization scale. Such a scale, distinct from 
Q and /i, may be useful when one or more invariants obey strong ordering. Here, however, we 
shall consider "fixed-angle" scattering configurations, in which the parameter Q sets the scale 
for all invariants, up to numbers of order unity. With this in mind, we will simplify Eq. fl2.4|l 
somewhat, and pick Q' = fi, that is, equal factorization and renormalization scales. Both the 
soft and jet functions then depend on ^^(yU^) only, and we will suppress their Q'^ dependence, 
now expressing the same amplitude as 

Ml^(^P.,^,a,{^i'),e^ = jm(a.(/.^),e)5S(A,^,a.(/i^),.) 

X iff (^/5„^,a,(/i2)^ , (2.6) 

that is, we suppress dependence on those variables that are set to unity by our choice of scales. 

Clearly, any jet-soft-hard factorization of the sort described above is unique only up to finite 
factors in the various functions. There is an additional ambiguity between the jet and soft 
functions at the level of a single infrared pole per loop in dimensional regularization. In the 
remainder of this section, we will provide specific definitions for the jet and soft functions that will 



enable us to define and resum them unambiguously, and which will be useful in our calculations 
below. We begin with the jet functions. 



A. The jet functions and the Sudakov form factor 



The factorization ()2.4|1 holds for any exclusive amplitude, including the elastic, or Sudakov, 
form factor. A very natural definition of the jet functions is, therefore, the square root of the form 
factor j^. Here, we will choose the case of the elastic scattering form factor with a color-singlet 
source, and spacelike momentum transfer. Reverting to the general case of jet momentum scale 
Q'^, not necessarily equal to the renormalization scale, this is 
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Q 



12 



(2.7) 



Below, we shall take /i as the MS renormalization scale, /i^ = /ig exp[— £:(7£; — ln(47r))]. With this 
choice, we may rely on the explicit form of the quark spacelike electromagnetic Sudakov form 
factor in D = A — 2e dimensions. A similar definition may be given for gluon jets in terms of 
matrix elements of conserved, singlet operators. In either case, the all-orders expression for the 
(square root of the) resummed form factor, organi zing all pol e terms, and implicitly specifying 



all finite terms of the jet defined as in Eq. (j^ZZj), is 3, |24| 



Q_ 
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exp 




1 

+ 2 



/i2 



f2.8) 



where we use a notation for the running coupling that emphasizes its re-expansion in terms of 
the coupling at fixed scale /x. For our purposes below, we shall need only the "leading" form of 
the running coupling. 



a. 



) 

^ / n=0 



(2.9) 
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with the one-loop coefficient 

/?o = yC^-^T^^F. (2.10) 

In the expression for the jet functions above, the choice Q'^ = /i^ can be imposed trivially. The 
functions /C^, and 7^ are anomalous dimensions that can be determined by comparison to 
fixed-order calculations of the Sudakov form factors for quarks and gluons. These form factors 



are now known in QCD up to three loops 
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26, 



27[. Notice that the coupling in the argument 



of /Cl^l is fixed at fi, so that the integral of this term alone is not well-defined at = even for 
D ^ 4. This apparent divergence, however, is cancelled by contributions from the upper limit of 
the /i^ integral of the anomalous dimension 7]^, and relates the latter to /C'*' order-by-order in 
perturbation theory. We will provide explicit expansions for these functions in Appendix 1X1 



B. The soft function 

We will broadly follow Ref. ^ in the definition of the soft function for partonic amplitudes, 
although we will modify certain details in the construction. The fundamental observation of 
Ref. i is ..t tKe soft funcUon, describing co.o. e.cKa,,e be.wee. the jets, is independent of 
coUinear dynamics, and may be constructed from an eikonal amplitude, that is, the vacuum ex- 
pectation of products of ordered exponentials. For each external parton of flavor fi, we introduce 
a nonabelian path-ordered phase operator, 



$[{^l((T',a) = Pexp 



(2.11) 



where f f ~ /3f is a four- velocity. For specific calculations at two loops, it will be useful to choose 
these velocities to be slightly timelike, 

0<t;2<l. (2.12) 

The "opposite moving" velocity projects out the large component of ff . The gauge field A^-^^'^ 
is a matrix in the representation of parton i. In the construction of the soft function, we will 

eventually take all vf 0, or equivalently, ff (3^. In perturbation theory, the operators 

\f-] \f-] 

$1// (cx3,0) and $1/- (0, —00) respectively generate outgoing and incoming eikonal lines in the Vi- 
directions. The eikonal sources couple to gluons at vertices in the color representation of parton i. 



An essential feature of these diagrams is that they are invariant under rescahngs of the velocities, 
Vi avi. 

We are now ready to construct eikonal multi-point amplitudes from products of ordered expo- 
nentials, tied together by the same color tensors, cl that appear in the expansion of the partonic 
amplitudes, Eq. ()2.4p . For the 2^2 case, 1 + 2 — 3 + 4, this gives 





X -oo).„.,$[^l(0, -oo),„,, |0) . (2.13) 



Such a product is gauge invariant. The eikonal amplitude, or web function, W depends in general 
on both the invariants Vi-Vj and the invariant lengths vf. The basic observation of Ref. ^ is that 
all potentially collinear divergent ratios factorize from dependence on wide-angle radiation for 
eikonal as well as partonic amplitudes. We can use this factorization to isolate the soft function 
systematically, using only calculations in the eikonal approximation. 

Because of the factorization of collinear singularities, such dependence is universal, depending 
only on the number and flavors of the external jets. In particular, as observed above, form 
factors, with two external lines and trivial color flow, generate the same collinear dependence. 
Thus, all collinear dependence cancels in the ratio of our four-point eikonal amplitude Wi and 
the product of two eikonal form factors, just as in the ratio of the four-point partonic amplitudes 
to the corresponding form factors. We shall deflne Sli by this ratio. Notice that information on 
color flow is not affected at all by the eikonal jet functions, which like partonic jets, are diagonal 
in color. Thus, we define 



where as above the velocities j3i are the lightlike limits of the Vi. The denominators are eikonal 
versions of the elastic form factors, defined with incoming velocities Vi and outgoing Wj, where 
vf = vf = v^ and Vi ■ Vi = uq, with uq a constant of order unity, independent of i, namely 





(2.14) 




(2.15) 
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This form factor generates the square of the colhnear poles associated with the eikonal jet of 
flavor i in Wu, and hence the soft function ()2.14|1 is free of collinear divergences. We may thus 
take the lightlike limit for the velocities to define the soft function in the ratio. 

Equation ()2.14|1 allows us to compute the soft function, once we determine how to choose the 
variable Uq, so that we may match the eikonal calculation to the partonic amplitude. We can 
determine the correct choice as follows. 

We first re-express Eq. fl2.6|) for the partonic amplitude, converting it into an expression for 
the soft function as a ratio analogous to Eq. ()2.14|) . 



This simple result enables us to set the scale Uq in the definition of the eikonal form factors of 
Eq. ()2.14|1 . In Eq. ()2.1(ij) . S can depend on the velocities only through the ratios (3i ■ (3jQ'^/ 
When S is calculated in this way from the ratio of partonic quantities, Q sets the scale of all 
momenta in the amplitude, and /i, the factorization scale in Eq. (j2.6|) . may be reinterpreted as 
the momentum transfer in the form factors that define the jet functions. When calculated from 
the eikonal ratio, on the other hand, S depends only on the variables f3i ■ Pj/uq. To match the soft 
function computed in eikonal approximation with the partonic amplitude, we need only require 

This relation will be used in our explicit calculations later. We are now ready to provide an 
all-orders expression for the soft function, analogous to Eq. ()2.8p for the jet functions. 



C. Resumming the soft function 

We will use the MS scheme for renormalization throughout. Before renormalization, all of 
the purely eikonal amplitudes discussed in the previous subsection give (only) scaleless integrals 
in perturbation theory. Such integrals vanish identically in dimensional regularization. In fact, 
these functions are only nontrivial because of renormalization, with every infrared pole resulting 
from the subtraction of a corresponding ultraviolet pole. This is the case whether or not W is 
collinear-regulated by introducing masses for its eikonal phases. 
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Thus, for both the web function W and the soft function S, we have (suppressing indices) 



SZe = 1 = Zs,iasifi),E)Sll, (2.18) 

and similarly for the eikonal form factors in the ratio ()2.14|1 . Both S and W are therefore defined 
entirely by their anomalous dimension matrices, 



1\ d{ZA)KJ 



= (^.-),,/3(.,.)^^, (2.19) 

which are given in any minimal scheme by the residues Z^^\, with A = Wf or S'f, of single 
ultraviolet poles in 1/e, at fcth order in the expansion 

CO , OO h ^ / T \ 

2^ = 1 + E (t) = E (t) E (I) ■ P-^") 

k=l k=l n=l ^ ^ 

Then, for example, from the one-loop bare integrals we find the one-loop anomalous dimension 
from the residues of the one-loop single ultraviolet poles, 

r« = -2Z«, (2.21) 

where F^'' is the nth order coefficient of (as/vr)" in F^. Similarly, to order 0{al), after one-loop 
renormalization we find the two-loop anomalous dimensions from the two-loop single poles, 

T^J^ = -AZ^^\. (2.22) 

From the definition of 5*, Eq. (j2.14j) . the soft anomalous dimension matrix is found from the 
matrix for the corresponding eikonal amplitude by simply subtracting the anomalous dimensions 
for the eikonal jets. We denote the latter by F2^(Mo/f ^, as), and write 



i^Sfjj = lim 



Wf,iJ\ — 5— '"'^ ) -^iJ Z^^2 ( 



(2.23) 



In Tsfi all sensitivity to collinear dynamics, and therefore to the choice of f^, is cancelled, and 
the coefficients depend only on the invariants (3i ■ j3j. 
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The matrix renormalization group equation for the eikonal amphtude Sjj^ is then 

'ft • AO' 



- r 



Si,IJ 



a. S 



[f] 

JK ' 



from which we can solve directly for 5 as a path-ordered exponential, 



Mo 



P exp 



1 



no 



(2.24) 



(2.25) 



where boldface (with a subscript for flavor flow) indicates a matrix. In summary, the matrix of 
anomalous dimensions, and hence the soft matrix itself, can be computed order-by-order purely 
from eikonal diagrams. 



III. THE JET FUNCTIONS TO TWO LOOPS 

In this section, we expand the jet functions in the factorized amplitude ()2.6p to fixed (second) 
order in in a form that is convenient for comparison to explicit partonic calculations. 

To determine the jet anomalous dimensions, as well as to use the resummed forms of the jet 
and soft functions with fixed-order calculations, we re-expand the running coupling in terms of a 
coupling at fixed scale. It is important to do so consistently in dimensional regularization, using 
the explicit form for the running coupling, Eq. ()2.9|) . It will also be convenient to use Eq. ()2.8|) 
as a starting-point to isolate the truly universal pole terms in the logarithm of the jet function, 
separating them from the finite terms. To this end, we introduce the notation 

n=l ^ m=l n=l ^ ^ 

= E^^{as{fi'),e) + e^Has{fi'),e), (3.1) 

in terms of the coupling ^^(/i^) at fixed scale fi. As in Eq. (j2.6|) . we set the jet factorization scale 
Q' = fi. The pure pole terms in Eq. p.ip have been expanded at each order as 

n+l p[j]{n)/ x 

^,Kn)(,),^^!!L_(£), (3.2) 

m=l 

while the functions e'*^ '■"-'(e) absorb all terms that remain finite for e = 0, order-by-order in as- 
The coefficients -Em and the functions e^'^^'^^^e) are determined, of course, by the expansions 
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of the functions 7^-, K, and which depend in general on the definition ()2.7|) of the jet. This 
separation, however, ehminates the remaining arbitrariness in choosing the form factor by defining 
a "minimal" jet, consisting of the exponential of pole terms only. 



In this notation, we rewrite our basic factorization, Eq. ()2.(j|l . in "minimal" form as 



(3.3) 



2 ' 



(3.4) 



where we have absorbed the (color-diagonal) finite factors into the perturbative definition of the 
short-distance function 



"Uf (^A, as{iJ^),e^ = exp 



Hf^ (A,^,«s(/i^)) . 



(3.5) 



We will also find it useful to write this expression in the color state notation of Eq. (|2.2p . as 



,as{fx^),e \Hi) 



(3.6) 



where again the matrix structure of the soft function is denoted by boldface and where we treat 
in the notation of Eq. (Q. 
Before this refactorization, the logarithm of the full jet function at two loops is given by 



lnjM(a>^)..) = i -(^)(^7r" + V"'>( 



+ 1- 

TT 



8 [ab^^ +^ * 



To determine the coefficients Em in the minimal two-loop jet function, we only need to expand 
the functions ("■)(£), 



2 + £ 



(3.7) 



^WW(£) = ^MW + e^WW'(0) + 



(3.^ 



12 



Explicit forms for these anomalous dimensions can be found in Appendix El In terms of these 
quantities, we readily find that the full single-pole terms in the logarithm of the jet function are 
given at one and two loops by 

1 



E. 



'i\ (1) 



El 



i\ (1) 



E. 



E. 



«](2) _ 3/5o [j](i) 



128 



Tk 



[i] (2) 



El 



i] (2) 



[<?] (2) 



32 

i-ic(2) + C(3) 



— C^Cf 
lo 



+ CpTpUF 

lo 



[9] (2) 



32 ^ 



346 11 
■^ + yC2 + C3 



+ —CaTftif 
16 



961 
216 



64 
27 



yC(2) - yC(3)' 



3 



C pTpHF ' 

16 



(3.9) 



where for we give the explicit expressions for the quark and gluon cases. Notice that the 

full single pole term includes a contribution from the running of the finite term at one loop, 
which appears as an 0{e) contribution in 



IV. EIKONAL AMPLITUDES AT ONE AND TWO LOOPS 

We begin this section with a calculation of the soft anomalous dimension matrix for quark- 
antiquark elastic scattering at one and two loops, in terms of a specific color basis j3|l, and then 

nn 

discuss the representation of the matrix in the color generator notation of [13|, [IJ] . We will see 
that the basic result of our calculation, the proportionality of the one- and two-loop matrices, 
applies to a much wider class of processes. 
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(a) 



(b) (c) 



FIG. 1: One-loop diagrams that contribute to T^g^ 
A. 2^2 eikonal diagrams at one loop 

Here we will present the calculation for one loop corrections to W, Eq. ()2.13|1 . for quark- 
antiquark scattering, and by using Eq. ()2.21|1 we will derive the corresponding one-loop soft 
anomalous dimension matrix. Representative one-loop diagrams are shown in Fig. ^ One can 
write the amplitude for any diagram D as 

MD = FnxCDj, (4.1) 

where is the corresponding Feynman integral in dimensional regularization and Co,i is the 
color tensor. We will refer to Fd as the velocity factor below, because it absorbs all dependence 
on kinematic variables. To uniquely define the normalizations of the velocity factors, and hence 
the color tensors, we define them to equal the corresponding integrals for the scattering of eikonal 
lines that couple to the exchanged gluons via color-independent "abelian" vertices. In particular, 
we absorb into the velocity factors the (—1) associated with a gluon coupled to an eikonal line in 
the anti-quark representation. Note that this separation of color and velocity factors is possible 



even if the eikonal lines are in the adjoint representation, 
comparison to results expressed in the formalism of Ref. 



'his method will facilitate our eventual 



3- 



Consider the left-hand diagram in Fig. ^ which we will call a "t-channel diagram" , referring to 



the pair of eikonal lines to which the gluon is connected. We will follow Refs. [28|, |29[ , and express 
the velocity factor as an integral in configuration space. For an arbitrary one-gluon correction to 
a phase operator of the form of Eq. ()2.13|) . such a correction is given by 

Ft = iigfin' [ dx, [ dy^D^^ix-y), (4.2) 

JCi JCj 

where integration is performed over the positions of gluons on the paths of the Wilson lines, Ci 
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and Cj. For the lines in Eq. ()2.13|) these paths are specified by 



Ci = Vif3, Cj = Vja, (4.3) 

where a and /3 run from — oo to (0 to oo) for an incoming (outgoing) path. For the t-channel 
diagram shown in Fig. Ufa), for example, where t = {pi — ps)^ = {p2 — P^)^, we may have 
{^,j} = {l,3} or {2,4}. 

In Feynman gauge the coordinate-space gluon propagator, in dimensional regularization with 
D = 4 - 2e, is given by 3 

= ^^2-. (a;2 • (4-4) 

Using this expression in Eq. (j4.2j) . we have 

/■oo fO 

.e\2 



Ft = {igp.") / da dl3v'^Df^^{vja — Vil3)Vj 

Jo J~oo 

= i^9l^'nv..v,)^—^ da dp (4.5) 
47r Jo Jo [{Vja + ViP)^ - te\ 

As observed above, all such integrals vanish in dimensional regularization, since they are scaleless. 
The contribution of each such velocity-dependent integral is given by its counterterm, equal to 
its infrared pole and hence to the negative of its ultraviolet (UV) pole. Of course. Ft may be 
evaluated as a momentum-space integral with equivalent results. 

In order to isolate the (single) UV pole in Eq. ()4.5|) . we apply an infrared cut-off for the 
integral by introducing a small parameter A with units of mass. This can be effected simply by 
inserting ^^(1/A — a) in Eq. ()4.5j) . The a and P integrals are then easily related to a single integral 
in terms of 2; = ^.^"^^ (see Eq. ()B3|) in Appendix |B)) . We find 



F^ = -(-)(^Yni-e)- fdz -. (4.6) 

The single UV pole term in this expression is given by j3[ 

F'-^\vi, Vj) = - (^^) ^-fij coth , (4.7) 
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where 



cosh 7,,- = ^^i;^ . (4.8) 



vfv] 



Because there is only a single, overall infrared divergence in Ft, any such cut-off will give the 
same ultraviolet pole. 

In the high-energy limit {'jij ^ 1), we have 

i^r(^.^.) = -(f)^7..- (4.9) 

For {i,j} = {1,3} and {i,j} = {2,4} the answers are identical, in this 2 — > 2 process. In the 
high-energy limit we define 

7l3 = 724 = 2" , 7l4 = 723 = U , 7i2 = 734 = 'S' , (4-10) 

where 

T^lnT^Uln^'"^-"^ 



y2 I \ y2 



C/ = ln('^Uln^'"^"^-^ 



y2 J y y2 

with = v"^ for all i. The velocity factors for u- and s-channel diagrams are found by taking 
into account the extra minus sign associated with coupling to an eikonal line in the antiquark 
representation, as well as that from crossing substitutions, which change the sign of coth 7^^ from 
unity to —1 in the high-energy limit, 

Fu{vi,Vj) = - Ft{vi,Vj) , 

Fs{vi,vj) = Ft{vi,-Vj). (4.12) 

Here, {1,3} = {1, 4} and {2, 3} for the xi-channel diagrams, and = {1, 2} and {3, 4} for the 
s-channel diagrams. The function Fg has the same overall sign as Ft because it differs both by 
an antiquark connection and by crossing, while F„ has the opposite sign. 
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ci = \ / C2 ; 




FIG. 2: Color basis {01,02} for four-quark process 
In summary, the single poles for the velocity factors for the diagrams in Fig. ^ are given by 

Fr^-{^)1-S. (4.13) 

To construct the counterterms, of course, we must also compute the corresponding color tensors 
for each diagram. 

We will use C|"^ to denote the color tensor for the t-channel diagram in Fig. Q with color 
tensor q, i = 1, 2 at short distances. For the latter we choose the basis tensors shown in Fig. |21 
The coefficients of the color tensors absorb all overall factors not included in the velocity factors 
of Eq. Km . 

One can calculate these color tensors from the basic identity for the generators of SU(A'"c), 

^ ^ )r2ri {T )r'3r4 = o '^»'2r4'^»'3^'l ~ r) AT ^r'2ri^r3r4 • (4-14) 



2 '"^ 2N, 

a 



In the color basis given in Fig. El the color tensors of the t-channel diagrams are given by 

4*' = ^|]^^2 = ^^^2. (4.15) 

We will employ a similar notation below for other one-loop and for two-loop diagrams. Color 
tensors for the u and s channel diagrams are computed in a similar way with the results 

C^-r^-^f^, (4.16) 
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and 



2 



2N, 



■C2 • 



(4.17) 



We summarize these relations in matrix form by 



(4.18) 



J=l,2 



where the matrix element dj] specifies the mixing from color tensor c/ to tensor cj by the 
exchange of a gluon in channel a = t,u, s. An important identity that we will use below is 



(4.19) 



which we easily verify from the relations above. Note that this equality holds in an arbitrary 
representation. ^ 

The contribution of each diagram to the matrix counterterm is now found by the product 
of the corresponding color factor times the pole part of the velocity factors, 



a=s,t,u ' 



(4.20) 



in terms of the single-pole terms of Eq. (|4.13p and the color factors read off from Eqs. (j4.15|) - 



(|4.17|) . Given the counterterm matrix, we can evaluate r[^| by using Eq. (j2.2H) with the result 



^ ^iU-T) + 2CpS 



.(1) 

Wt 



\ 



\ 



(S-U) 

iT-U) ±.{U-S) + 2CpT ) 



(4.21) 



Exactly the same calculation gives Eg , the anomalous dimension for the eikonal jet function, 
defined as the square root of the eikonal singlet form factor, Eq. (j2.15p . In Eq. (j4.9|) . we simply 



^ Equation H4.19|l is equivalent in this case to the weU-known identity = in the color generator notation 

that we will review below. 
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let cosh ■jij — > uq/v^ = ji^/iQ'^v'^) (using Eq. ()2.17|) ). in the hmit v"^ 0. The one-loop result 



for parton i is then given by 



[i](l) po 



^ = :r a In 



2,, 2 / ' 



(4.22) 



with Ci = Cp for quarks and Ca for gluons. By using this expression and the definition for F^j, 
Eq. (12121, we find 



^ St 



(S-U) 

{T-U) ^^(U-S) + 2CfT J 



(4.23) 



where 



(4.24) 



After performing the subtraction of the jet functions, we set Vi — > Pi, and then use Eqs. ()2.3|) and 
()2.17|l to recast the result in terms of the usual Mandelstam variables, t, u and s. We notice that, 
as anticipated, all collinear logarithms, and hence sensitivity to our choice of collinear regulation, 
are absent in the soft anomalous dimension matrix, T^^^"*. 



B. 2 — > 2 eikonal diagrams at two loops 

Figure El shows the classes of topologically inequivalent diagrams that contribute to F^^ , when 
combined with their one-loop counterterm diagrams. One obtains the full set from all different 
combinations of external legs with these topologies. It is easy to see that the number of graphs 
for each inequivalent set is A^^ = 6, A'f, = 6, = 6, A'd = 12, = 12, Nf = 12, Ng = 4, 
Nfi = 24 and Aj = 3, which in total gives 85 two-loop diagrams. As in the one-loop case, we find 
anomalous dimensions from the combinations of velocity factors and color tensors. 

Consider first diagram {i), which is the only two- loop topology involving all four eikonal lines. 
Diagram (i) does not have a surviving single UV pole when we add its one-loop counterterms. 

Regarding the remaining cases, we consider first those diagrams involving two eikonal lines 
only, which we refer to as "2E" diagrams. Next, we will show that diagram (g) vanishes, which 
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(a) (b) (c) 




(d) (e) (f) 




(g) (h) 0) 

(2) 

FIG. 3: Two-loop diagrams that contribute to T^/^ 

we consider a very important result. Finally, we will calculate all the contributions from the 
surviving "3E" diagram type {h). In this case, we will find that the diagrams, although non- 
vanishing, reduce to the product of one-loop diagrams, and thus do not contribute to the two-loop 
anomalous dimensions. 



1. The 2E diagrams and T\ 

In the 2E diagrams, (a), (6), (c), (rf), (e) and (/), the gluons connect to only two of the 
four eikonal lines in W . These same diagrams also contribute to the two-loop cusp anomalous 
dimension, Fg' ^'^\ and their single UV poles are well known Q]. We review their velocity factors 
here, because they are needed for the two-loop anomalous dimension matrix. Additional details 
are given in Appendix O 

The color factors of the 2E diagrams are proportional to the color-mixing matrix elements for 
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single t-channel gluon exchange, (ij^, defined in Eq. ()4.18|) . This is manifestly the case for the 
individual diagrams (c), (d), (e) and (/). For the sum of diagrams (a) and (6), it relies on the 
result 0, that the single-pole terms in the velocity factors of these two diagrams are the 
negatives of each other. The net color factor for the (a) and (6) single poles is then proportional 
to the commutator of two generators, which allows it to be expressed in terms of the one-loop 
color factor, as Ca d}j\- We can thus present the contributions of all the 2E diagrams in terms of 
the Sj], with a = s,t, u. 

In terms of the factors the two-loop counterterms^ for the diagrams (a), (6), (c) and (/) 
in the high-energy limit are, analogously to the one- loop velocity factors, Eq. ()4.13j) . 

2 C, 



JI 



■T^ C(2)^_C(3)' 
6 2 2 



7{C) 



JI 



''J I 



2 2e V 36 



3V 5 



and 



JI 



a 



> ,[t] Ca 1 



9 



- + (i-C(2))r 



2^3 

"6 



+ 



T2 

— -T+a2) 



(4.25) 



(4.26) 



with T the logarithm of 2vi-v^/v'^, as in Eq. 1)4.111) . In Eq. ()4.26|) . the second term in square brack- 
ets gives the result of those numerator terms that are proportional to before the integration. 
(See Appendix O) The entire T-dependence of these terms cancels against the contributions 
from diagrams [d) and (e), which are also proportional to w| before integration and are given 
individually by 

2 



7{d) 



7{e) 



JI 



JI 



2 2 



^) d 

TT 



[t] 
JI 



C(2)' 



(4.27) 



^ These results, of course, require that we combine these diagrams with the corresponding one-loop counterterms 
for their divergent subdiagrams. Notice that diagram (6) in Fig. O does not have a one-loop counterterm since 
it does not have a divergent subdiagram. 
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The combined t-channel contribution from the six diagrams (a), (6), (c), (d), (e) and (/) to the 
soft anomalous dimension matrix is found by adding Eqs. ()4.25|1 - ()4.27|1 ^. 



ylt]] —of yia+b) I 7(c) I r, 



I 7(e) I yif) 



JI 



2e Ktc 



''JI 



V36 2 J 9 

+ ^ (C(2)-C(3)) 



In 



2vi ■ V3 



Tx ) \ 2 (a./vr) 4^ 



(C(2)-C(3)) , 



(4.28) 



where the second hne recalls a standard notation [2d| for the quantity K 
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(4.29) 



The result ()4.28j] includes a factor of two for the other t-channel exchange, between lines 2 and 
4. 

Analogous considerations, of course, apply to diagrams with pairs of s- and w-channel 2E 
diagrams. Together with the t-channel diagrams, they contribute to the two-loop anomalous 
dimension matrix for W according to Eq. ()2.22|) . 

-2eFr- 

2 ^ ' ~ 



r(2i?)(2) 



JI 



i=s,t,u 



+ 5j,CAa(C(2)-C(3)) 



K 
~2 



rii; + 5j/C^a(C(2)-C(3)) , 



(4.30) 



where we have used the identity ()4.19|1 , and where r[^| is the same one-loop anomalous dimension 
given in Eq. ()4.21|) . 

In a precisely similar manner we find for the two-loop form factor (cusp) anomalous dimension 
for partonic representation 



9i 
4 



K In 



+ Ca(C(2)-C(3)) 



(4.31) 



Note that one needs to multiply Eqs. I|4.26|l and 14.27|l by 2 because for these diagrams there are two ways of 
attaching the gluons to the eikonal lines. 
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As at one loop, we combine Eqs. ()4.3U|) and ()4.31|) in Eq. ()2.23|) . in order to find the contribution 
of tlie 2E diagrams to tlie two-loop soft anomalous dimensions for scattering. 

It is now clear that the two-loop soft anomalous dimension matrix inherits from the 2E dia- 
grams a factor of K times the one-loop anomalous dimension matrix. The result is, 

rr^'^ = f r«, (4.32) 

with rlj^ the same one- loop anomalous dimension given in Eq. (j4.23|) . All velocity-independent 
terms in the Z^^"^ that are not in F^^^'' cancel in Eq. ()2.23|) against the corresponding finite terms 
from the eikonal form factors in the two-loop soft anomalous dimension, along with all coUinear- 
singular dependence. This is important, because the constant terms depend in general on the 
eikonal approximation and our choice of collinear regularization. At the same time, we have now 
used a// the collinear-singular dependence in the Sudakov anomalous dimensions, Eq. (j4.3H) . to 
cancel the In f ^-dependence of the 2E diagrams of W . The 3E diagrams, represented by {g) and 
{h) in Fig. El have no remaining subtractions. The combination of these classes of diagrams must 
therefore be free of collinear singularities at the two-loop level. 



2. Vanishing of three- gluon diagram with three eikonal lines 

Now let's show that diagram (g) in Fig. El vanishes. Up to overall factors which play no role, 
the velocity Feynman integral for a generic three-gluon diagram can be written as 

11 1 11 



F{va,vb,vc) = J (fki(fk2 



1 

{ki + k2y + ie 



vb ■ ki + ie va ■ k2 + ie vc ■ {ki + ^2) + ie kf + ie k^ + ie 
va -vbVc- {ki - /C2) +va-vcVb- {ki + 2/^2) +vb -vcVa- {-2ki - /C2) 

(4.33) 



where the term in square brackets is the three-gluon vertex momentum factor. Here va, vb and 
Vc are three different eikonal velocities. We take lightlike v\ = v"^ = 0. We can then expand any 
momentum p'^ as 

= —VB-P + ^VA-P + P^, (4.34) 

Va -Vb Va- vb 
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with the transverse components, satisfying va ■ px = Vb ■ Pt = 0- 
For use in the integral we introduce the variables: 



. va-vc , 
4i = vb ■ ki, 



vb ■ Vc , 



(4.35) 



va - Vb va- Vb 

We introduce these variables into the integral by using va and to define light-cone coordinates, 

1 



dkjdk- 



■ d{vB ■ ki) d{vA ■ ki 
Va ■ Vb 

Va ■ Vb 



(va ■ vc){vb ■ Vc) 



diidrji 



(4.36) 



so that 



2 _ o (^A • ki){yB ■ kj) _ ^2 ^ 2 '"a-vb 

va-vb ''^ {va- vc)ivB ■ Vc) 



(4.37) 



When we change variables in F to the .^'s and r^'s, we find 

2 



F{va,Vb,Vc) 



Va ■ Vb 



{va ■ vc){vb ■ Vc) 



X- 



I Xldiidr^id^'-^^T 

i=i 



' (vA-Vc){l'B-Vc) 



' {VA-VC){VB-Vc) 
1 1 



(^1 + 6) (^1 + V2) - iki,T + k2,T) + ie 
1 



6 + ier]2 + ie ^1 + 6 + ^71 + ^72 - vc,t ■ {ki + ^2)^ + 

6 - 6 + ^1 - ^72 - Vc,T ■ {ki - k2)rp + ^1 + 2^2 - 2771 - r]2 



0. 



(4.38) 



The integral vanishes because the numerator is antisymmetric under ^1 772, ^2 ^ Vi 
ki,T ^ ^2,r; while the product of the denominators is symmetric. Notice that group factors play 
no part in this argument. This result is therefore very general and applies to any 2 n process 
with lightlike velocities. 



3. Exponentiation of the remaining 3E diagrams 



The only remaining class of diagrams is illustrated by diagram {h) in Fig. 01 Along with its 
companions found by permuting the eikonal lines, we can refer to these as "3E" diagrams, since 
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(a) (b) (c) 



FIG. 4: a-c) Pairs of 3E diagrams. 

they are the only nonvanishing diagrams with gluons connected to three eikonal hnes. There are 
24 such 3E diagrams, 8 of them with s and t channel gluon exchanges, 8 of them with s and u 
channel gluons, and finally 8 of them with t and u channel gluons. They come in pairs as shown 
in Fig. m 

We now show that the analysis of the previous subsection regarding the three-gluon diagrams 
leads to a very interesting result for the remaining 3E diagrams as well. In this case, the dia- 
grams do not vanish, but reduce to products of one-loop diagrams^. They therefore provide no 
contribution to the two-loop anomalous dimension matrix. 

Each such 3E diagram contains one eikonal line with two gluons attached to it, which we 
label as vc- The two lines having one gluon attached are labelled va and vb- We will consider 
a pair of 3E diagrams that are related simply by exchanging the order in which the two gluons 
attach to vc, as in Fig. Ufa) for example. The two diagrams have differing color and momentum 
structures, but we can rewrite their sum as the sum of one term with symmetric color and 
momentum integrals, plus a second term with antisjTiimetric color and momentum integrals. In 
the spirit of the discussion above, we suppress the color matrices held in common and write 

^ah{vA, Vb, Vc) = J^j^ff'^\vA, Vb, Vc) + J^ir^^^^V^, Vb, vc) , (4.39) 

^ The reduction of a different class of multi-loop eikonal diagrams, namely the 2E diagrams of ladder type, to 
powers of one-loop diagrams, was previously observed in Ref. [sj to hold to all loop orders. 
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where the subscripts a and h refer to the color generators on the fc-eikonal, contracted with 
generators on the va or vb hne, respectively. Consider first the antisymmetric term, which is 
given by 



Vb ■ ki + it Va ■ k2 + ie kf + ie fc^ + 



vc ■ ki + ie vc ■ k2 + ie 



(4.40) 



Vc ■ {ki + k2) + ie 

(In the color-generator notation described in section IIVDI below, the color operator associated 
with this antisymmetric term takes the form [T^ ■ Tc , Tc ■ T^a\-) '^^^ same change of variables, 
Eq. ()4.H5j] . leads to an expression that is again manifestly antisymmetric under the relabelling 

6,2 ^ ^72,1, ki^T ^ k2,T, 

:F^r''''^\vA,VB,vc) = \ inT^-T^n) 

1 r ^ 1 

{va ■ vc) {vb ■ vc) J l\ ^ {v^-va){vl.va) ^^^^ - Kt + 

111 1 

X 

6 +teri2 + le ii + rji- vc,t ■ fci,T + ^e ^2 + V2 - vc,t ■ k2,T + «e 
+ vc,T ■ {k2 - ki)j, 



^i+^2 + Vi+V2- vc,T ■ {ki + k2)j, + ie 
= 0. (4.41) 

The entire color-antisymmetric part of the infrared region thus vanishes whenever the eikonal 
approximation is valid, and the cancellation is exact for the eikonal amplitudes we consider here. 

Turning to the symmetric term, we need only use the eikonal identity l/[x{x + y)] + l/[y{x + 
y)] = l/{xy) to rewrite it as the product of the two lowest-order single-gluon exchange diagrams, 

T^r\vA.VB.vc) = \ {nTa + T,n) 

X i d^'k ^ ^ ^ 



Vb ■ ki + ie Vc ■ ki + ie ki + ie 

r — ■ r — -72^- (4-42) 

Va - k2 + le Vc ■ k2 + le fc| + le 
(In the color-generator notation described in section IIVDI the color operator associated with 
this symmetric term takes the form {T^ ■ Tc , Tc ■ T^i}.) These diagrams have the correct 
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color and kinematic structure to represent the two-loop terms in the exponentiation of the color 
matrix of one-loop infrared poles. They are precisely cancelled in the two-loop soft function by 
the corresponding products of one-loop counterterms. 

We conclude that the entire two-loop anomalous dimension is due to the 2E diagrams and is 
given by Eq. ^.3^ , in which we may remove the superscript {2E), to obtain 

rg' = frw. (4.43) 

We have thus determined that the two-loop anomalous dimension color-mixing matrix is related 
to the one-loop matrix by the same factor that relates the one- and two-loop Sudakov anomalous 
dimensions, A{as)- Evidently, the next-to-next-to-leading poles in amplitudes with color ex- 
change are generated by the same exponentiation of "webs" as for the elastic form factor 
Additionally, we note that in the "bremsstrahlung" or CMW scheme ^4 \ , this contribution, along 
with the corresponding term in the cusp anomalous dimension, is absorbed into a redefinition of 
the strong coupling, which effectively boosts the strength of parton showering. 

C. Expansion of the soft function 

To relate the soft anomalous dimension to fixed-order calculations, we expand the resummed 
soft function, given as a path-ordered exponential in Eq. ()2.25p . to order 0{a'j.). The result is 



n'+iBM-rS'. (4-44) 



where in the second equality we have used Eq. ()4.43p for the two-loop anomalous dimension 
matrix. Combining this result with the second-order minimal jet function, Eq. ()3.7|) . in the 
formula for the factorized amplitude, Eq. ()3.4|) . we will derive a result to compare directly with 
the pole structure of explicit two-loop calculations. 
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D. Iterative color-matrix form 



Given the one- and two-loop anomalous dimension soft matrices (I4.23|l and ()4.43|1 . and the 
expansion of the quark jet function, as in Eq. ()3.7|1 . we can use the factorized amplitude, Eq. ()3.4|1 . 
to calculate all infrared and collinear poles at order for quark-antiquark scattering. The result 
will be a set of coefficients of the specific basis tensors in color space that we have chosen. Fig. |21 
In this basis, we can perform threshold resummation for jet and other cross sections. 

To compare to explicit calculations at the two-loop level, however, and to generalize to higher 
numbers of external partons, it is convenient to make contact with a somewhat different notation, 
in which the color interactions of soft gluons is represented by a color matrix for the insertion 
of a gluon on external line i, with a generator in the color representation of that parton i, 
whose color-matrix (rather than generator) indices are summed against those of the lower-order 
amplitude that is "dressed" by this soft gluon. In the notation of Eq. ()2.2j) above for the color 
content of an amplitude, the action of the generators may be made explicit as the action of a 
vector, with indices in the adjoint representation; for example for i = 1, 

\mM ^M^^ ff3^^9l^as{i^'),e) 5. (T'^),^,^ , (4.45) 

-Ia!i,r2... \ H J 

where & = ±1 absorbs minus signs associated with antiparticles and crossing. In the convention 
of Ref . , 5j = 1 when % is an eikonal line representing an outgoing quark or gluon, or incoming 
antiquark; 5j = — 1 for an incoming quark or gluon, or outgoing antiquark. Defined in this 
way, the vector color generator matrices obey the fundamental relation Tj = 0, which is an 
expression of gauge invariance. The TjS are conventionally normalized to Tj ■ Tj = Cj, with 

A gluon exchanged between two parton lines i and j produces the product Tj ■ Tj, which acts 
on an amplitude in a fashion precisely similar to Eq. ()4.45|) . This notation allows for a convenient 
iterative expression for color exchange due to soft gluon exchange, without requiring an explicit 
choice for the color basis. 

The color-space notation above may be applied to the computation of the soft function as 
well as to the amplitude itself. Consider the soft function at a single loop, determined by the 
one-loop soft anomalous dimension. As we have seen, the latter is built up from the contributions 
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of soft gluon exchanges between pairs of eikonal lines. From each exchange, the contribution to 
the anomalous dimension is found from Eq. ()2.21|1 . where the one- loop single-pole term in 
equals the one-loop UV pole term computed from the corresponding diagram. 

Each diagrammatic contribution, then, is proportional to a product Tj ■ Tj acting on the 
lower-order amplitude, multiplied by the result of the eikonal integral. Referring to Fig. ^ the 
relevant single-pole coefficients are given in Eq. ()4.13|) . The action of T^^^ on the color tensor is 
the sum of all such terms, with a subtraction for the jet anomalous dimensions; this subtraction 
is proportional to the identity matrix in color space, as in Eq. ()2.23j) . This gives 



where Sij = {pi + pj)"^, with all momenta defined to fiow into (or out of) the amplitude. For the 
four-parton case above, S12 = s, S13 = t, and so forth. The overall 1/2 compensates for double 
counting in the sum. To derive the final result, we have used the explicit forms of the 5jS described 
above, as well as the identities J2i Tj = and Tj ■ Tj = Cj (in the quark-scattering case, all 
Ci = Cp)- In this notation the color identities enforce the cancellation of the collinear-sensitive 
ln(l/f^) terms. 

Identical considerations apply to the two-loop case. The nonvanishing anomalous dimension 
matrix is again a sum of diagrammatic contributions, corresponding to gluon exchange processes 
involving two eikonal lines only. As we have seen, these contributions have the same color- 
generator structure, Tj ■ Tj, found at one loop. The 3E diagrams have a more complicated color 
structure, but they do not contribute to the two-loop soft anomalous dimension matrix. 

To be more specific, we saw that diagram {h) in Fig. |21 can be organized into antisymmetric 
and symmetric color structures, which can be represented as commutators and anti-commutators 
of one-loop color structures, of the form [Tj ■ Tj , T^ ■ T^] and {Tj ■ Tj , Tj ■ T^}. Note that the 
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\M,) 



\M,) 



(4.46) 
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antisymmetric quantity can be written as, 

[Ti ■ Tj , Tj ■ Tk] = V [T", Tj] = -i TJT^ . (4.47) 

As the final form sliows, it is totally antisymmetric under permutations of the three eikonal lines. 
This is also the form of the color factor for the other type of 3E diagram, the three-gluon diagram 
{g) in Fig. El 

As emphasized above, the velocity factors multiplying both commutator and anti-commutator 
structures vanish. (In the case of the anti-commutator, the vanishing occurs after adding the 
one-loop counterterms.) Nevertheless, we display the commutator in Eq. ()4.47j) . because it has 
occurred in the literature before. We will encounter such terms below in our analysis of explicit 
two-loop calculations, and show how they are consistent with the specific solution for the soft 
anomalous dimension, Eq. ()4.44|) . in which this combination of color generators does not appear. 

E. Generalizations 

The analysis given above applies far beyond the 2^2 quark-antiquark scattering amplitude. 
When the soft anomalous dimension is expressed in terms of color generators, as in Eq. ()4.4(j|l 
at one loop, and using this equation and Eq. ()4.43p to do so at two loops, the result is slightly 
less explicit than, say, Eq. ()4.23|) . but it is much more general. When we generalize from quark 
and antiquark to gluon lines, and when we add more partons in the final state, the only change 
in our considerations above is to change the color generators Tj, and sum over more variables i 
and j. The eikonal momentum integrals that give rise to the coefficients of the generators are 
the same for any choice of parton pairs or triplets. 

In these terms, the two-loop results organized in Eqs. ()4.46|) and ()4.43|) are not limited to 
quark-antiquark scattering, but apply to the scattering of any flavor combination. Furthermore, 
these relations are by no means limited to 2 2 scattering, and apply to any 2 ^ n process. 



as in multi-jet production. These resu 
resummations in hadronic scattering 



ts, therefore, are a step toward threshold and related 
at the level of next-to-next-to-leading logarithm. 
At present, however, for the purposes of resummation we must still rely upon the explicit 
form of the matrix as in Eq. ()2.25|) to generate the amplitude at arbitrary orders. Anticipating 
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further applications, it will be useful to investigate flexible choices of color basis, perhaps based 
on the trace notation described, for example, in Ref. (this point was noted in Ref. We 
reserve these considerations for future work. 



V. SINGLE POLES AT NNLO 

In this section, we combine the expansions of the jet and soft functions in the "minimal" 
factorized amplitude, Eqs. ()H.4j) and and give an explicit expression for infrared poles to 

two loops, including single poles. We go on to compare these "postdictions" of the two-loop 
single pole terms to the results of explicit calculations, and verify that they agree. Traditionally, 
these results have been presented in a form proposed some time before by Catani and we will 
briefly review this formalism and relate it to the two-loop expansion of our resummed expressions. 

A. Two-loop poles from the factorized amplitude 

Here, as above, we adopt the notation f{as) = Xlnl'-'^s/^)"/''"^- ^^^^ notation, we can 
express the Born and one-loop amplitudes for process f in terms of the factorized jet, soft and 
hard functions of Eq. as 




where for the jet functions we have used the minimal form (j3.3|) . In the second equality for 
|A4f^''), we have used explicit expressions for the jet functions and the soft matrix, the latter 
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from Eq. ()4.44|) . Using these results, we find for the two- loop amplitude 
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Now both the and the s["^ are given by sums of pure poles in e. As a result, their squares 

and products all begin at At two loops, then, the single-pole terms that multiply the 

Born amplitude |A^[°^) in Eq. ()5.3|) are given entirely by the single poles in E'W^^) and S[^\ From 
Eq. ()3.3|) for the jets, and Eq. ()4.44|) for the soft matrix, these poles are found from the coefficients 
of the soft and jet anomalous dimensions, 
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where we have separated the double- and higher-order pole terms from single-pole terms that 
multiply the Born amplitude, followed by single poles times the one-loop hard scattering, and 
finally the two-loop hard scattering. 

From Eq. ()3.9|) . the two-loop single-pole terms that multiply the Born amplitude l^Af^) in 
Eq. ()5.4p are given by 
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where we recall the notation of Eq. (j3.8|) for the coefficients ^^("^(e:). Given that the one- and 
two- loop have been known for a long time, and that we have just calculated the two- 

loop soft anomalous dimension matrix, this expression provides an explicit form for the intrinsic 
two-loop single poles in dimensionally regulated amplitudes. 
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Note that a redefinition of Sj to include a non-pole term would both change the definition of 
\'Hf^) at one loop in Eq. ()5.2j) . and introduce single-pole terms into the (E'W'^^^ + Sj^'')^ contribution 
to the two- loop expression Eq. ()5.3|) . As we shall see below, the Born-times-single-pole terms 
remain invariant under this shift only if the shift commutes with the F^^'*. We therefore need all 
of the expansion ()5.4|1 . in order to make contact with the results of explicit two- loop calculations 
at the single-pole level. 



B. One- and two-loop amplitudes in Catani's notation 

To compare to existing calculations, we now review the notation of Ref. Q], in which they are 
normally presented. We first observe that in this notation, amplitudes are organized in powers 
of (as/27r), rather than (as/vr). We will distinguish this trivial difference below by a prime in 
the color states, as |A^[(")) = 2"|A<f(")). 

In this formalism, the single- and double-pole structure of one-loop amplitudes is expressed 
in terms of the color generator operators introduced above. 



-Si 



(5.6) 



with ^^(l) = — 7_B the logarithmic derivative of the Gamma function. Here 7^ = (3q/2 and 
= Ca for gluons {i = g), and 7j/Tf = 3/2 for z = g or g. The poles of the one- loop amplitude 
in color-state notation are then represented as 

l^l^W) = \M'/o)) + . (5.7) 

The explicit relation to the resummation formalism at one loop is found by expanding in 
powers of e, 
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Taking into account the overall factor of two from the expansion in asjlT^-, the pole terms in 
Eq. ()5.7|1 are thus identical to those in Eq. ()5.2|1 . The matrix ll^"* generates as well explicit 
/^-dependent finite contributions contained in If^^'^'^, which in the "minimal" factorization are 
absorbed into the one-loop hard function |7if). The one-loop infrared finite amplitudes are 
related by 

= + i^^w) . (5.9) 

This is an example of a finite shift of the sort mentioned above, which redefines the finite function 
at one loop. 

At two loops, Ref. [3l predicted the fourth- through second order poles in terms of the 
generators I^^'', and absorbed the then-unknown single pole contributions in terms of a color 
operator H^^-*, 



^"^r^r^ ^\^-) + Hfl.) . (5.10) 



6^^^(1)1(1 - 2e) //?o 
T{l-e) \2e 

The two-loop amplitude is then organized as 

= |A^;(0)) +l(l)(e) I^JW) + |^;(2)fin) . (5.11) 

In the intervening years, the color generators Hj have been determined by matching to the 
single-pole structure of explicit two-loop QCD scatte ring amplitude calculations, for example 
99 ~^ 99) ~^ 99) ~^ QQi and e^e" — > qqg |l^ [ll[ [iJQ, Is^- Here we follow Ref. j^l and 
write 

Hf)(.) = l|5:/7f + Hf^)| + 0(.), (5.12) 

where we split the single-pole factor into a color-diagonal term, which can be represented as a 
sum of constants Hi for each external parton i, and a matrix HJ; that includes all color mixing. 



This matrix can be written as 
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Hf ) = I Yl fa.a^as V' Tf T«3 In ( In lu ( ) , (5.13) 

(i,j,k) ^ ^^'^ ^ \ Skt / \ Sij 
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where the sum is over distinguishable but unordered triplets of external lines k). We note the 
similarity to the color structure from the 3E diagrams, Eq. ()4.47jl . We emphasize that this form 
has been obtained directly only for processes with at most four partonic legs. In Ref. [2^ it was 
also shown to be consistent with the proper coUinear behavior of the 2 — n gluon amplitudes. 

It is these expressions that we will compare to the two-loop expansion of the factorized am- 
plitude, Eqs. ()5.4p and ()5.5p . Rather than provide explicit expressions at this point for the 
constants H^'^'^ from [2^, we will derive below expressions relating the constants H^'^'^ to the iet 
anomalous dimensions (in MS scheme). Here we will find useful an identity found in Ref. (8fil |. 
The matrix Hf ^ in Eq. will emerge from our results for the two-loop soft anomalous di- 

mension matrix, plus the effects of a one-loop finite color-mixing term. We now turn to this 
exercise. 



C. H^^) from the anomalous dimensions 

Inserting the definition oi I\ , Eq. ()5.10|) . into Eq. 1)5.111) and expanding to the accuracy of 
we readily find 

= i (i;"(e))'|M;<°'> 
I (l!"(2.) -!["(.)) + I<"(2.) + Hf + W»'> 
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We will relate this expression to the single-pole result from the factorized amplitude, Eqs. ()5.4|) 
and 

The single-pole terms in Eq. ()5.14j) that multiply the Born amplitude come from two sources: 
the (If^^)^ operator on the first line, and the terms in the square brackets on the second line, 
which also include finite corrections indicated by +0{6^). 

To make contact with the expansion of the resummed amplitude, Eq. ()5.4|) . we first separate 



35 



the poles in each of the terms, according to Eq. ()5.8|) . 
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where in ij^-*^" we isolate the finite terms from I^^'* that multiply the Born amplitude. Comparison 
with Eq. ()5.4p requires further that we commute the soft anomalous dimension matrices with 
poles to the left of the finite amplitudes, and that we also re-express 

|_A/f/(i)fin^ in terms of 

using Eq. ()5.9p . The first step, in particular, leads to an additional commutator contribution at 
the level of the single poles times the Born amplitude, 
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Here we have organized the expression just as in Eq. ()5.4|1 . starting with the square of one-loop 
pole terms, two-loop second- and third-order poles, and then first-order poles, all times the Born 
amplitude, followed by poles times the one-loop hard amplitude and the finally the two-loop hard 
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part, 

We are now ready to compare this expression to the two- loop single pole terms of Eq. ()5.5p . 
Higher-order poles can easily be checked in a similar manner j9|. 

Consider first the matrix parts of Eqs. ()5.5|1 and Recalling the factor of 4 associated 

with changing from the coefficient of (as/vr)^ to (as/2vr)^, we see that the K F^^^'' term is identical 
in the two expressions. Consistency then requires the remarkable result that the commutator of 
I^^^*^"^ in Eq. ()5.8|) with the one-loop soft anomalous dimension in Eq. ()4.46|) precisely cancel the 
two- loop hI^-* function as defined in Eq. ()5.13|) . 

l?'\T^^=-\lif\ (5.18) 
In fact, a compact calculation, given in Appendix El shows that Eq. (j5.18p indeed holds for 



the explicit matrix Hj given in Eq. ()5.13j) . for arbitrary 2^2 processes, and also for 2 
processes where all particles are identical. For those processes with five or more partons for 
which the quantities 7i/T^ are not all identical, the commutator is more complicated, as can be 
seen by inspecting Eq. ()5.8|) . and as discussed in Appendix |D| Because we know the anomalous 
dimension matrix for all these processes, however, Eq. ()5.18|) can be turned around and taken as 
a definition of the corresponding matrices Hj; . Recently, the soft anomalous dimension matrix 
r^^^'' for 2 — s> 2 processes was computed Q| by making use ofjust this connection to H^^"*, along 
with the explicit results for h[^'' for quark-quark scattering lfll |. 

The remaining, color- diagonal, single-pole terms in Eq. ()5.16p are found using the values of 
the one-loop quantities ii^W(^) given in Eq. dSH), and the form of H^^ in given in Eq. ()5.12|) . 
Then the single-poles times Born amplitudes of Eq. (j5.16|) are given by 
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where we have suppressed dependence that contributes only at the level e^. The comparison of 
Eq. ()5.19|) with the expansion from the resummed amplitude, Eq. ()5.5|) . is now trivial. We simply 
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appeal to the striking identity noted explicitly by Ravindran, Smith and van Neerven 
in our notation is written as 



which 



4 32 ' 4 

(2) 



where is given in Eq. ()3.9|) . In Ref. [36| this expression was observed to imply a close 

(2) 

relationship between the constants and the form factors. We now see that, aside from color 
mixing, all the single-pole terms are identical to those in the form factors. Indeed, the precise 

(2) 

terms relating the to the single-pole residues of the elastic form factor are present simply 
to cancel a set of "extra" single-pole terms generated from the expansion of I^^^ in the two- 
loop amplitude. As in the case of the color-mixing anomalous dimensions, we can also consider 
Eq. (Ih:^ definition of the constants H^'^K 

In summary, we have shown that the full single-pole structure of the two-loop amplitudes can 
be reconstructed from the same anomalous dimensions that determine the next-to-next-to-leading 
poles of the factorized jet and soft functions at all orders in pertubation theory. This relation, 
and the explicit forms of the anomalous dimensions, hold for partonic scattering amplitudes with 
arbitrary numbers of external lines. 



VI. CONCLUSIONS 

We have extended the factorization and resummation formalisms for exclusive amplitudes in 
QCD to next-to- next-to-leading poles in these amplitudes. The same anomalous dimension ma- 
trices, calculated here directly for the first time at two loops, control a variety of resummed cross 
sections at NNLL. These calculations generalize the determination of the Sudakov anomalous 
dimensions to nontrivial color mixing. 

We verified the formalism and anomalous dimensions by showing that they allow us to repro- 
duce the very nontrivial color and momentum structure of single infrared poles at next-to-next- 
to-leading order for 2^2 processes in the literature. 



The calculation of the NNLO soft anomalous dimensions opens the door to 
mation at next-to-next-to-leading logarithm for multijet cross sections 



threshold resum- 
7|. Perhaps our 



38 



most striking result is the discovery that the two-loop soft anomalous dimension matrix is ob- 
tained from the one-loop matrix simply by multiplying by Kas/ {^ti), where K is the constant 
given in Eq. ()4.29|1 . This is exactly the same property obeyed by the scalar Sudakov or "cusp" 
anomalous dimension. Aside from its intrinsic interest, this relation will make possible next- 
to-next-to-leading logarithmic resummation formulas in a closed form, since it will be possible 



to diagonalize the two-loop anomalous dimension matrix independently o: 



coupling 



33|, I34I, using the same color eigenvectors found at one loop |4 l6l ItI. ISTL l38|. 



the run ning of the 



Our analysis applies not only to 2 ^ n processes relevant to hadronic colliders. In addition, 
it applies to the inelastic scattering of a parton by a color-singlet source, as in deep-inelastic 
scattering, and to the creation of arbitrary numbers of partons from a color-singlet source in 
leptonic annihilation. It will clearly also be of interest to extend this analysis to massive external 
lines. 
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APPENDIX A: ANOMALOUS DIMENSIONS 

In this appendix, we provide the low-order anomalous dimensions entering the jet function, 
as defined in Eqs. ()2.7|1 and ()3.7j) . We give the nth-order coefficients 71^^"^ /cW(") ^^d in 
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an expansion in powers of Q;s(/i^)/7r, 
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The resuhs for ^H(") -^ere obtained from Ref. [23], which also contains results through three 
loops. We shift the gluonic expressions by terms proportional to /5-function coefficients, which 
take into account the effects of the renormalizing the operator C^^^G""^" , as explained in Ref. 3^. 
Because we only quote results through two-loop order, some of the results for could also 



have been extracted from the two-loop quark electromagnetic form factor 
Higgs boson amplitude 



25[ and from the gg 
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APPENDIX B: ONE-LOOP VELOCITY FACTORS 



1. Basic integrals 



Consider the one- loop t-channel diagram shown in Fig. ^a). The velocity factor is given by 



Ft = (igfi^)'^ (f 1 ■ V3) / da dj3 D{v3a + vi(3) 
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We will use the following change of variables 



Jo [(t^3« + vi(3Y] 



a + P = T] , a = ZT] , 

with Jacobian rj. For IR regularization we impose a < j 
< z < 1. In terms of the new variables we have 



Ft = {ig^i 



.2 1 



(Bl) 



(B2) 



V < j^- ^ Also note that 
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with z' = ^ — 1- From the expansion of the above expression, the single-pole term and the finite 
part of the one-loop diagram are given by 

2 ' 



In ( ^ ) + ln(7re 
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where we have defined the following integrals 



h{vi,V3) = dz 



dz' = , 

{vs z + vi zf Jo {vs + vi z'y ' 

, ln(t;3 + vi z' f 
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with z = 1 — z. 



^ For these one-loop diagrams there is one overall IR divergence since all the collinear singularities factorize. 
Therefore it is sufficient to restrict only one of the giuon attachments. 



41 



2. Evaluation of Ii and 

We are evaluating eikonal diagrams derived from external Wilson lines. By looking at the 
usual momentum-space expressions for the amplitudes one can easily see that all these diagrams 
are scale independent in the eikonal velocities Vi. With this property in mind we can simplify 
the evaluations of the integrals by choosing vf = 1 without loss of generality. 

In order to evaluate Ji(fi,f3) we use the following change of variable j28| 

^2, ^ V^^+^MEflL , (B6) 



which gives 
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dz 2 {v^z + vizy 

From this change of variable it is very easy to see that 



Therefore we get 
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It is easy to verify that the first and last terms in the right-hand side of the final expression 
cancel, for example by changing variables to u = y + e~^^^ in the first term and u = y + e"^^^ in 
the last. This leaves us with 



ImiVl,V3) 
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In the high-energy limit 7^1, where 3> e~^, one easily finds 
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Following the same steps one also gets 
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APPENDIX C: VELOCITY FACTORS FOR 2E DIAGRAMS 



We begin our analysis of the 2E diagrams with the diagram that has a three-gluon vertex, 
Fig-Et/)- We follow Refs. 28|,|39| and write the three-gluon vertex as 

Vp(A;, /, -k-l) = V^,,{k, I) + D^,p{k, I) , (CI) 

where 

V^up{k, I) = {21 + k^^Qyp + 2kpg^y — ^kj^g^p , 



Dp,up{k, I) — luOfip + ^)p9p.v 



(C2) 



and where k and / are the loop momenta. Indices u and p attach to the ^3 line. The diagram 
resulting from Vp^up is proportional to f| before integration. We also note that the contributions 
of diagrams of Fig. Efrf) and Fig. |3fe) are entirely proportional to f | before integration, since a 
single gluon propagator attaches twice to the same eikonal line. These ^3 contributions turn out 
to cancel each other in the high-energy limit. We give the result for the Vp,yp contribution below. 
The contribution resulting from the D^^p piece for the diagram of Fig. Ef/) is given by, 
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where we have suppressed factors of ie in the denominators. One can evaluate the above ex- 
pression in either momentum or configuration space. We will not review the derivation of the 
following result 



e)'^4^{o.3Coth,.(-l+^C(2))4«0.3)} 
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^0 sinh 7i3 - sinh \ sinh ip 

We analyze I2 in order to get the high-energy behavior of this amplitude. We start by writing 
I2 as 

I2 = Icth-l + II 5 
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and where 
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Note that Icth-i is exponentially suppressed in ip when ip ~ 713. However, for small ip the factor 
"'"^^^^^ 2 + C(e-^i3). Therefore we can rewrite hth-i as 
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We evaluate ki, k2 and separately, starting with 
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Combining these results one finds 

Wi = C(2)7i3. (C13) 

Now let's look at the remaining integral in Eq. ()C6|) . Ji. After some trivial algebra one can 
rewrite Ji as 

h = k4 + k5, (C14) 
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Finally is given by 
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by using same kind of manipulations. Combining the above results one finds 
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Using Eqs. ()C13|) and ()C18|) . one finds for the asymptotic behavior of I2 

By using this result in the expression for the single-pole term, W2E 3g-D' 
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The contribution of the V^^p piece to the diagram in Fig. Of/) is given by 
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One can analyze the high-energy asymptotics of I3 in a similar way as above, with the 
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Combining Eqs. (IHM . KJ2T\i and and letting 713 = T, we find 
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which is the result given in Eq. ()4.26p . 

The amphtudes for diagrams (d) and (e) of Fig. IHl are given in Ref. (2^J and the high-energy 
asymptotics is obtained with a similar analysis. The results are given in Eq. ()4.27|1 . As mentioned 
above, these contributions cancel the V^i^p contribution from diagram (/), which is enclosed by 
the second set of brackets in Eq. ()C24|) . 

Finally, let us look at the crossed-ladder diagram in Fig. ^b) . The velocity factor in configu- 
ration space is given by 
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It is not difficult to show that the single-pole part of the crossed-ladder velocity factor is precisely 
the negative of that for the uncrossed-ladder diagram in Fig. IHl^a). Therefore in the combination 
of the two diagrams the single-pole part of Eq. ()(^25j) is multiplied by the difference of the 
respective color factors. Although the individual color factors are not proportional to the one- 
loop factor d^jj, their difference evaluates to <ij/C*A/2- The following result for the combination 
of diagrams (a) and (6) can also be found in Ref. j28|,^ 

W^clV,* = - (v)' 4^ Ye ^0*^^13/3^,^3) , (C26) 

where we define 

^713 

/a iv,,V3)= di/ji/j (713 - ^) coth ^p. (C27) 
^0 

One can investigate the asymptotic behavior of J3 in a way similar to that presented for the 
3-gluon vertex diagram in Fig. Et/). One obtains the result 

h{v„vs) = ^ + ^713 - ^ + 0(e-^-) . (C28) 
2 2 



By using the above relation in Eq. (jC26jl we find 

Letting 713 = T, this is the result given in Eq. ()4.25|1 . along with the result for diagram (c 

^ Needless to say we can evaluate the integrals in momentum space and get the same result. 
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APPENDIX D: THE COMMUTATOR OF I 



(l)fin 
f 



AND r 



(1) 
Sf 



The task of this appendix is to evaluate the commutator 

(1) 



(l)fin ^(1) 
f ' 5f 



appearing on the left- 



hand side of Eq. ()5.18p . Note that the pole parts of l\ can be identified with F^^ , via Eq. ()5.8|) . 
Writing out the 0{e^) parts of I^^'*'^'^ with nontrivial color structure, the commutator of the finite 
and pole parts of becomes, 



(l)fin ^.(1) 



' ^ St 



k l^k 



-In^ 



-Ski 



Ik , 

^ In 

^k 



Ski 



■ T.) In 



-Si 



(Dl) 



where in the second equality we introduce notation to separate the terms with three logarithms 
(Cs^f) from those with two (C2,f). 

In the case where all external lines are gluons, or all are quarks and/or antiquarks, all the 
ratios 7fc/T^ in the left-hand side of the commutator are equal. This term is then proportional 
to r^^,\ and C2,f vanishes. This argument does not apply, of course, to mixed processes, such as 
qq gg. For the latter case, however, and for any other 2^2 process, we may use the color 
conservation identity Ylk Tfc = and the simplicity of the kinematics to show that C2,f vanishes. 
The argument is simple, and may be given for the case qq gg without loss of generality. In 
this case, we may take = 1, 2 in Eq. ()D1|) to correspond to the incoming quark and antiquark, 
and we consider just these terms in the double-logarithmic part of the commutator in Eq. ()D1|) . 



We focus first on the terms with prefactor 'Jq/Tg = 3/2. (The same argument applies to the 
remaining terms, with prefactor jg/Tg.) These terms are proportional to 



V V T, ■ T, In 



k=l l^k 



-Ski 



i=l j^i 



(D2) 



This commutator would vanish if the sum over index k were extended to k 
can by done by observing that J2k Tfc = implies that, for example. 



3 and 4. But this 



Ti ■ T2 = T3 • T4 + ^ (Tl + - Tl - Tl) 



(D3) 
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where the squared terms commute with all combinations of generators. At the same time, we 
have Si2 = S34. As a result, in the left-hand term of the commutator (jD2jl . we may make the 
replacement 



Analogous reasoning for each of the terms in the sum over k and / in Eq. ()D2|) shows that for 
2-^2 scattering the sum of the missing terms with A; = 3, 4 is identical in the commutator to 
the sum from k = 1,2. Inserting the missing terms, at the price of an overall factor of 1/2, the 
two entries of the commutator become identical and it vanishes. This argument, of course, is 
heavily dependent on the specifics of 2 ^ 2 scattering. We know of no general argument that 
would eliminate all double-logarithmic terms in the commutator in 2 ^ n processes; indeed such 
terms are generically present. 

We now consider the triple- logarithmic terms in Eq. ()D1|) . 



where in the first equality we have introduced the notation bki = ln^(/i^/(— s^i)) = hk and 
ttij = ln(/i^/(— Sjj)) = aji. In the second equality in Eq. ()D5|) we have identified the nonvanishing 
terms in the commutator, for which one and only one pair of generators is matched between the 
two entries of the commutator. Because the scalar products are symmetric, there are four ways 
in which this matching may occur, for fixed indices i ^ j ^ k. Finally, the third equality shows 
the result of performing the commutator explicitly for the generators on the j line. This form is 
reminiscent of the color structure of 'tlf'\ Eq. ()5.13|) . although the triple-logarithmic momentum 
factors are different, and depend on the renormalization scale, /i. 

To make contact between Eq. ()D5|1 and the explicit expression ()5.13|) for tlf'^ , we convert the 
sum of unequal choices of i, j and k into a sum over distinguishable triplets, denoted {i,j,k). 
For each such choice, there are six permutations of the indices i, j and k in the final expression 




(D4) 





(D5) 
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of Eq. ()D5|) . These can be thought of as three cychc permutations, which leave the structure 
constants the same, but change the momentum factors, and three more (exchanges of i and k for 
fixed j, plus cyclic permutations), which change the sign of the structure constants, and change 
the kinematic factors. 

Following this path, we define 

'^[kjji] = ^kj — bij 0.jk (D6) 

and rewrite Cs^f as 

C3,f = 2 ^ faia2a:i T"^ T"^ C[kjji] + C[ik,kj] + C[ji,ik] ■ (D7) 

A straightforward calculation shows that all of the /i-dependence cancels in this expression. 
Relabelling the indices, and using the antisymmetry of the structure constants, we derive 

C3,f = E /™TrTf T- In(^) In(^) In (^) 

{i,j,k) ^ ■''^^ \ -^ki/ \ ^ij / 

= -\ Hf ^ , (D8) 

which establishes the result of Eq. (j5.18|) . We emphasize that, unlike our demonstration that C2,f 
vanishes, this result holds for an arbitrary 2 — > n process. 
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